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Introduction: Figure 1 shows a radar image and 
topography for flat-topped, steep-sided “pancake” 
domes on Venus. At least 145 such domes have been 
identified on Venus [1] and are thought to be volcanic 
in origin [2], Based on analysis of the dome surfaces, 
[3] suggested that only the late stage surface fractures 
are preserved, indicating entrainment and annealing of 
fractures during emplacement, consistent with a basal- 
tic composition. 




Regio. (b) Magellan altimetry of a typical dome in 
the Rusalka region. 

Previous studies have attempted to model the em- 
placement of these domes [e.g., 4, 5] as a laminar axi- 
symmetric viscous gravity current of a Newtonian 
fluid [6, 7], There are, however, several significant 
fundamental shortcomings in the model presented in 
[6, 7], These are semi-analytic methods based on simi- 
larity solutions that satisfy unlikely boundary condi- 
tions (i.e., increasing volume flow rates). Even for a 
constant effusion rate, the suspect treatment of singu- 
larities leads to a flow front advance that goes as t 18 
[6,7] instead of t l/2 as shown earlier in [8] with singu- 
larities neutralized. Moreover, most formulations con- 
sider the fluid to be isothermal. 

Due to a host of physical processes, such as the 
cooling, formation, fracturing, and entrainment of a 
cmst, the bulk rheology of the domes during emplace- 


ment is unclear. Although a Newtonian flow rate is 
often assumed for simplicity, other rheologies are also 
admissible, particularly Bingham (e.g., [9]) and empri- 
cally derived flow rates (e.g., [10, 11, 12, 13]). Baloga 
et al. [10] indicate that letting n = 1 in (1), below, re- 
sults in a flow that glides on a thin basal layer, essen- 
tially a limiting case of the Bingham flow rate. Typi- 
cal values of n for empirically derived flow rates for 
terrestrial mass flows cluster around 1.5. 

Here we present a method for inferring bulk rheol- 
ogy from topographic profiles of domes. A widely 
referenced (but little known to the volcanology com- 
munity) contribution by Babu and Van Genuchten [8] 
presents an innovative perturbation solution to the ra- 
dial formulation of the so-called Boussinesq equation 
for fluid flow. This approach allows accurate and 
nonsingular approximations to be made near the origin 
to avoid this problematic singularity and allows for 
solutions to within 0.1% of the flow front. The most 
significant advantage of the Babu and Van Genuchten 
approach is the ability to distinguish a flow’s rheology 
based solely on the shape of the flow surface. 

Approach: The equation describing the radial ex- 
pansion of a fluid with a gently sloping free upper sur- 
face is given by the continuity expression 


\ d ( n dh\ dh 

\rh m — =— (1) 


r dr \ dr J dt 

where r and h are the radial and vertical variables, re- 
spectively, t is time, and the value of n is indicative of 
the fluid rheology. The radial variable, r, can have 
values ranging from 0 (origin) to ry(t) at the flow front. 

To solve the boundary value problem when the 
steady state volume flux, Q, at the origin is given, 
equation (1) is transformed to an ordinary differential 
equation. To find the similarity solution by perturba- 
tion, the parameter, e, is defined as 


e = 1 t(n + l) 
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A new independent variable, x, and dependent vari- 
able P, are defined as 
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The constants in (4) are found to be 0/ ~ 0.5775, and 




0 2 ~ 0.7455, by requiring that x — 1 at the the flow 
front, r f . Equation (1) can now be rewritten as 
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From [8], the similarity solution of (5) is 
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Application: The examples shown here use the 

approximation from [8] to evaluate P(x) for 0 < x < 
0.1. A Romberg numerical integration scheme was 
used to evaluate all the integrals in (6) and (7) to find 
values of P(x) for 0.1 < x < 1. For illustration pur- 
poses, a dimensionless volume flux of Q = 50 is used. 

Assuming a Newtonian rheology as in [10, 11] 
with n — 3 in (1), Figure 2 shows the solution for a 
radially spreading fluid at four times. The small peak 
near the origin is a numerical artifact of the approxi- 
mation used near the singularity. However, even with 
this small artifact, the overall “shape” of the flow sur- 
face, as well as the aspect ratio at the final time, is very 
similar to the example shown in Figure lb for a steep- 
sided pancake dome on Venus. 



Fig 2. Axially symmetric Newtonian fluid flow pro- 
files at dimensionles times, t = 5, 35, 70, 140 from 
left to right, respectively. 


Perhaps the most intriguing aspect of the solution 
described here is that the topographic profiles for dif- 
ferent rheologies have distinctly different shapes (Fig 
3). The most distinct differences are between the limit- 
ing Bingham “basal glide” and Newtonian flow rates. 
The basal glide flow has a concave upward surface 
that gently slopes from the origin to the flow front. 
This shape is not at all consistent with the altimetry 
data presently available for the pancake domes (e.g. 
Fig lb). However, the Newtonian flow rate produces a 
concave downward surface that is characterized by an 
essentially flat top and a steep flow front, very similar 
to that observed on Venus. 
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Fig 3. Upper free surface at t = 140 for three vol- 
ume flow rates as defined by the n values shown. 
All other boundary conditions are the same. 

Conclusions: This new approach has broad poten- 
tial for modeling lava domes with a simple algorithm 
that does not require complicated finite difference 
schemes. The approach has the added advantage that 
it can be used to distinguish between basic flow rates 
for fluid emplacement. Based on this very preliminary 
assessment, the pancake domes on Venus are consis- 
tent with a Newtonian bulk rheology. 
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